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ABSTRACT 

In this paper, we compare the observationally derived black hole mass function 
(BHMF) of luminous (> 10^^ - 10^^ erg/s) broad-hne quasars (BLQSOs) at 1 < 
z < 4.5 drawn from the Sloan Digital Sky Survey (SDSS) presented in Kelly et al. 
(2010), with models of merger driven BH growth in the context of standard hierar- 
chical structure formation models. In the models, we explore two distinct black hole 
seeding prescriptions at the highest redshifts: "light seeds" - remnants of Population 
III stars and "massive seeds" that form from the direct collapse of pre-galactic disks. 
The subsequent merger triggered mass build-up of the black hole population is tracked 
over cosmic time under the assumption of a fixed accretion rate as well as rates drawn 
from the distribution derived by Merloni & Heinz (2008). Four model snapshots at 
z = 1.25, z = 2, z = 3.25, z = 4.25 are compared to the SDSS derived BHMFs of 
BLQSOs. We find that the light seed models fall short of reproducing the observa- 
tionally derived mass function of BLQSOs at Mbh > W^Mq throughout the redshift 
range; the massive seed models with a fixed accretion rate of 0.3 Edd, or with accretion 
rates drawn from the Merloni & Heinz distribution provide the best fit to the current 
observational data at z > 2, although they overestimate the high-mass end of the 
mass function at lower redshifts. At low redshifts, a drastic drop in the accretion rate 
is observed and this is explained as arising due to the diminished gas supply avail- 
able due to consumption by star formation or changes in the geometry of the inner 
feeding regions. Therefore, the over-estimate at the high mass end of the black hole 
mass function for the massive seed models be easily be modified, as the accretion rate 
is likely significantly lower at these epochs than what we assume. For the Merloni & 
Heinz (2008) model, examining the Eddington ratio distributions fsdd, we find that 
they are almost uniformly sampled from fsdd = 10~^ — 1 at 2: ~ 1, while at high 
redshift, current observations suggest accretion rates close to Eddington, if not mildly 
super-Eddington, at least for these extremely luminous quasars. Our key findings are 
that the duty cycle of SMBHs powering BLQSOs increases with increasing redshift for 
all models and models with Pop HI remnants as black hole seeds are unable to fit the 
observationally derived BHMFs for BLQSOs, lending strong support for the massive 
seeding model. 



1 INTRODUCTION 



Demography of local galaxies suggests the that most galax- 
ies harbour a quiescent super-massive black hole (SMBH) 
in their nucleus at the present time and the mass of the 
hosted SMBH is correlated with properties of the host bulge. 
Observational evidence supports the existence of a strong 
connection between the growth of the central SMBH and 
stellar assembly of the host spheroid (Tremaine et al. 2002; 
Ferrarese & Merritt 2000; Gebhardt et al. 2000; Marconi & 
Hunt 2003; Haring & Rix 2004; Gultekin et al. 2009) and 
possibly the host halo in nearby galaxies (Ferrarese 2002). 
This strong correlation is suggestive of co-eval growth of the 
BH and the stellar component likely via regulation of the gas 



supply in galactic nuclei from the earliest times (Silk & Rees 
1998; Kauffmann & Haehnelt 2000; Fabian 2002; King 2003; 
Thompson, Quataert & Murray 2005; Natarajan & Treister 
2009; Treister et al. 2011). Meanwhile, our observational un- 
derstanding of accreting black holes and their properties at 
high redshifts z > 6 is growing (see Fan et al. 2010 for an 
overview). Black holes appear to be copiously accreting and 
in place before the universe was a billion years old as evi- 
denced by the recent detection of the quasar a,t z — 7.085 
in the UKIDSS survey reported by Mortlock et al. (2011). 
Our access to the earliest redshifts is of course restricted 
to the brightest candidates, optical, and deep X-ray surveys 
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are needed to uncover the potentially highly obscured pop- 
ulation at these epochs (Treister et al. 2011). 

Simultaneously our theoretical understanding of the as- 
sembly history of black holes over cosmic time has been 
rapidly improving (see recent reviews by Volonteri 2010 and 
Natarajan 2011). With the increasing availability of obser- 
vations from earlier and earlier epochs, models can be effec- 
tively tested at extremely high redshifts. Recently, a more 
complete census of the accreting black hole population at 
1 < 2 < 4.5 has become available enabling the estimate of 
the mass function (MF) of these black holes (Kelly et al. 
2010, KIO hereafter). In this paper, we explore a range of 
theoretically viable models to fit this additional and new 
data-set. 

Observations currently provide strong constraints lo- 
cally, and hints at high redshifts (Mortlock et al. 2011; Treis- 
ter et al. 2011). With this additional anchor during inter- 
mediate epochs (z ~ 1 — 4), discrimination between models 
finally starts becoming possible as we demonstrate here. In 
particular, we are able to distinguish between seed black 
hole formation scenarios. This paper is the third in a se- 
ries (Volonteri, Lodato &Natarajan 2008, and Volonteri & 
Natarajan 2009, VLN08 and VN09 respectively hereafter) 
where we compare the late time evolution of SMBH seeds, 
focusing on the difference between a population of more nu- 
merous, but smaller mass BHs as seeds, versus a population 
of rarer, but more massive initial seeds. The outline of this 
paper is as follows: in Section 1, we briefly review the new 
observations and BHMFs derived therefrom; Section 2 pro- 
vides the context of the theoretical black hole growth mod- 
els followed by a description of how BHMFs are derived in 
Section 3; and finally we present the results of comparing 
observations to the models and the derived conclusions. 



2 BHMFS FROM OBSERVATIONS AND 
MODELS 

KIO derive an estimate of the BHMF of BLQSOs correcting 
for incompleteness and statistical uncertainties from a sam- 
ple of 9886 quasars at 1 < 2 < 4.5 from the SDSS. They find 
'downsizing' of BHs in BLQSOs, i.e. the peak of the number 
density shifts to higher redshift with increasing black hole 
mass peaking at « ~ 2. They report that as a function of 
back hole mass and Eddington ratio, the SDSS at z > 1 
is highly incomplete at Mbh < 10^ M0 and L/L-Edd < 0.5. 
The lower limit on the lifetime of a single BLQSO phase 
was estimated to be > 150 ± 15 Myr, with a maximum 
black hole mass of ~ 3 x W^°Mq. KIO also find that the 
Eddington ratio distribution peaks at L/I/Edd ~ 0.05 with a 
dispersion of ~ 0.4 dex implying that most BLQSOs are 
radiating nowhere near or at the Eddington limit. From 
their estimated lifetime and Eddington ratio distributions 
they also infer that most massive black holes spend a signif- 
icant amount of time growing in an earlier obscured phase 
consistent with models of self-regulated growth. The recent 
claimed discovery of a population of heavily obscured, co- 
piously accreting black holes at 2; > 6 by Treister et al. 
(2011) provides an interesting glimpse perhaps of the prece- 
dents of the SDSS BLQSOs population at lower redshifts. 
Note that this claim of a detection of a copiously accreting 
and obscured population of black holes at 2 > 6 on X-ray 



image stacking by Treister et al. (2011) as well as its signifi- 
cance have been disputed by Cowie et al. (2011) and Willott 
(2011). The issue can only be settled conclusively with more 
high redshift data that will hopefully be forthcoming. As 
noted in our earlier work, in particularly, VN09, with more 
comprehensive observational redshift coverage, several key 
assumptions in the modeling approaches can be tested. 

In this paper we focus our modeling efforts on four red- 
shifts of interest (2 = 1.25; 2 = 2; 2 = 3.25; 2 = 4.25; 
to match the redshift bins in KIO). At each redshift our 
models provide us with a sample of all SMBHs present at 
that particular cosmic time, and of the SMBHs that are 
actively accreting. From the SMBH mass and its Edding- 
ton ratio, /Edd, we can derive their bolometric luminosity: 
log(Lboi/ergs-i = 38.11 + \og{Mw/MQ) + log(/Edd). We 
apply a bolometric correction of 4.3 (as in KIO) and se- 
lect only quasars that are more luminous than the mini- 
mum luminosity determined by the flux limit described in 
Richards et al. (2006) , the parent sample of KIO. The thresh- 
old bolometric luminosities corresponding to each bin are 
as follows. At 2 = 1.25, log(Z/boi/ergs~^)=44.7; at 2 = 2: 
log(Lboi/ergs"^)=45.2; at 2 = 3.25: log(Lboi/ergs"^) 45.5; 
at z — 4.25: log(I/boi/ergs~^)=45.8. Finally, we assume that 
the fraction of unobscured quasars is 20%, based on La 
Franca et al. (2005). Here we do not explicitly apply the 
evolutionary model of La Franca et al. (2005), where the 
fraction of obscured quasars depends on both redshift and 
luminosity, as the redshift range we are interested is beyond 
that explored by La Franca et al. (2005), but we note that 
when we apply their evolutionary model to the redshift range 
z — 1 — 3, adopting the bolometric corrections of Marconi 
et al. (2004), we obtain consistent results. Recently Fiore 
et al. (2012) have also derived and published BHMFs for 
the entire active SMBH population as opposed to just the 
BLQSOs as done by KIO. 



3 BLACK HOLE GROWTH MODELS 

To trace the assembly history of black holes with cosmic 
time, Monte-Carlo realizations are performed to derive the 
merger histories of dark matter haloes. We also trace the 
formation and growth of embedded black holes as a function 
of cosmic time as outlined below. 



3.1 Black hole seeds 

To track the mass assembly history of black holes in the 
universe, we need to start with seeds at high redshift. In the 
standard picture, the assumption is that the remnants of the 
massive first stars (Pop III stars) provide the earliest seeds in 
the range of 50— 100 Mq . We note here that whether the first 
stars were indeed this massive has been called to question 
from the latest round of recent higher resolution simulation 
results where fragmentation occurs ubiquitously (Turk et al. 
2009; Greff et al. 2011; Davis et al. 2011). An alternate model 
for the formation of massive seeds from the direct collapse of 
pre-galactic disks is presented in Lodato & Natarajan (2007; 
2006) . In these models, there is a limited mass range of halos 
with a further narrow range in spins that are able to form 
seeds. However, contrary to the Pop HI case, massive seeds 
with M « 10^ - 10^ Mq can form at high redshift (2 > 15), 
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when the intergal actic medium has not been significantly en- 
riched by metals (IKoushiappas et al. [20041 : iBegelman et alj 



riched by metals (IKoushia.ppas 
I2OO6I : iLodato fc Nataraianll200' 



20071'). More details of this 



seeding model can be found in lLodato fc Natarajarj (|2006l . 
[20071), wherein the development of non-axisymmetric spi- 
ral structures drives mass infall and accumulation in a pre- 
galactic disc with primordial composition. The central mass 
accumulation that provides an upper limit to the SMBH 
seed mass that can form is given by: 



Mbh = mdMh 



for 
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A < Aniax = 771d(5c/8(md/jd)(rvir/Tgas)^ 



(1) 



(2) 



and AIbh = otherwise. Here Amax is the maximum halo 
spin parameter for which the disc is gravitationally unsta- 
ble, irid is the gas fraction that participates in the infall and 
Qc is the Toomre parameter. The efficiency of SMBH for- 
mation is strongly dependent on the Toomre parameter Qc, 
which sets the frequency of formation, and consequently the 
number density of SMBH seeds. Guided by our earlier inves- 
tigation, we set Qc = 2 (the intermediate efficiency massive 
seed model) as described in VLN08. 

The efficiency of the seed assembly process ceases at 
large halo masses, where the disc undergoes fragmentation 
instead. This occurs when the virial temperature exceeds a 
critical value Tmax, given by: 

2/3 

(3) 
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TTld 1 -I- MBH/mdMhalo 



where ac ~ 0.06 is a dimensionless parameter measuring the 
critical gravitational tor que above which the disc fragments 
iJRice. Lodato fc Armita ge 2005). 

To summarize the seeding model, every dark matter 
halo is characterized by its mass M (or virial temperature 
Tvir) and by its spin parameter A. The gas has a tempera- 
ture Tgas = 5000K. If A < Amax (see eqn. 2) and Tvir < Tmax 
(eqn.[3]), then we assume that a seed BH of mass A/bh given 
by eqn. ^ forms in the center. The remaining relevant pa- 
rameters are rrid — jd = 0.05, Ofc = 0.06 and here we con- 
sider the Qc = 2 case. In the massive seed model, SMBHs 
form (i) only in haloes within a narrow range of virial tem- 
peratures (10'* K< Tvir < 1.4 X 10* K), hence, halo velocity 
dispersion {a ~ 15kms~*), and (ii) for a given virial tem- 
perature all seed masses below nid M modulo the spin pa- 
rameter of the halo are allowed (see eqns. 1 and 3). The seed 
MF peaks at lO'^M©, with a steep drop at 3 x lO'^M©. We 
refer the reader to Lodato & Natarajan (2007) and VLN08 
for a discussion of the mass function (and related plots). 
Here we stress that given points (i) and (ii) above, the ini- 
tial seeds do not satisfy the local A/bh — cr relation, in fact 
the seed masses are not correlated with a, rather they are 
correlated with the spin parameter of the dark matter halo. 

In this paper, we once again contrast this model with 
a popular scenario that advocates the first black holes are 
the remnants of zero metallicity stars. We assume that one 
Pop HI star forms i n metal-free halos with Tvir > 2000 K 
([Yoshida et al.ll2006l 'l. We assume a logarithmically flat ini- 
tial MF, dN/d\ogM =const, bet ween lOAfa and e OOAfs , 
where the upper limit comes from lOmukai fc Fallal ((2003), 



and suppose that seed black holes form when the progeni- 
tor star is in the mass range 40 - UOAf© or 260 - QQQMq 
([Fryer et al.ll200ll ). The remnant mass is taken to be one- half 
the mass of the star. 

Both scenarios for black hole seed formation rely on zero 
metallicity gas. We model the evolutio n of metallicity by 
the 'h i gh feedback, best g u ess' m odel of IScannapieco et al.l 
(200v|). IScannapieco et al.l ([2003) rnodel metal enrichment 
via pair-instability supernovae winds, by following the ex- 
pansion of spherical outflows into the Hubble flow. They 
compute the co- moving radius, at a given redshift, of an 
outflow from a population of supernovae that exploded at 
an earlier time. Using a modificatio n of th e Press-Schechter 
technique, IScannapieco fc Barkanal ([20021 ) compute the bi- 
variate MF of two halos of arbitrary mass and collapse red- 
shift, initially separated by a given co-moving distance. From 
this function they calculate the number density of super- 
novae host halos at a given co-moving distance from a 're- 
cipient' halo of a given mass Mh that form at a given redshift 
z. By integrating over this function, one can calculate the 
probability that a halo of mass Mh forms from metal-free 
gas at a redshift z. When a halo forms in our merger tree 
we calculate the probability that it is metal-free (hence, it 
can form Pop HI stars) and determine if this condition is 
satisfled. 

Every halo entering the merger tree is assigned a spin 
parameter drawn from the lognormal distribution in Aspin 
found in numerical simulations, with mean Aspin = 0.05 and 
standard deviation a\ = 0.5 (Davis & Natarajan 2009). We 
assume that the spin parameter of a halo is not modified by 
its merger history, as no consensus exists on this issue at the 
present time. 



3.2 Black hole growth 

We evolve the population of SMBH seeds according to 
simple models of self-regulation with the host. The main 
features of the models ha ve been discussed elsewhere 
([Volonteri fc Nataraianll2009l ). We summarize below all the 
relevant modeling assumptions. SMBHs in galaxies under- 
going a major merger (i.e., having a mass ratio > 1 : 10) 
accrete mass and become active. Each SMBH accretes an 
amount of mass, AAf = 9xlO'^(cr/2OOkms"*)*-^''Af0, where 
(7 is the velocity dispersion after the merger. This relation- 
ship scales with the A/"bh — o" relation, as it is seen today 
(|Gultekin et al.ll2009l ): 



A/b 



1.3 x 10' 



f - V 

V200kms-V 



Mrr 



(4) 



the normalization in AAf was chosen to take into account 
the contribution of SMBH-SMBH mergers, without exceed- 
ing the mass given by the Mbh — o" relation. 

We link the correlation between the black hole mass and 
the central stellar velocity dispersion of the host with the 
empirical correlation between the central stellar velocity dis- 
persion and the asymptotic circular velocity as ct = v^jspi) 
of galaxies (?). The latter is a measure of the total mass 
of the dark matter halo of the host galaxy. We calculate 
the circular velocity from the mass of the host halo and its 
redshift. 

The rate at which mass is accreted scales with the Ed- 
dington rate for the SMBH, and we set either a fixed Ed- 
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dington ratio of /Edd = 1 (for Pop III seeds), /Edd = 0.3 (for 
massive seeds), or an accretion rate d erived from the distri- 
bution derived by Merfoni fc Heind (|2008) (we apply this 
model to massive seeds only). The empirical distribution of 
Eddington ratios derived by Merloni & Heinz (2008, MHOS 
thereafter) is fit by a function in log(I/boi/iEdd)- The fit- 
ting function of the Eddington ratio distribution as a func- 
tion of SMBH mass and redshift, is computed in 10 red- 
shift intervals (from 2 = to z = 5) for 4 different mass 
bins (6 < log(MBH/M0) < 7; 7 < log(MBH/MQ) < 8; 
8 < log(AfBH/A'/Q) < 9; 9 < log(MBH/M0) < 10), and then 
fit with an analytic function which is the sum of a Schechter 
function and a log-normal (A. Merloni, private communica- 
tion). The Eddington ratio distributions are then normal- 
ized to unity at every given mass and redshift. We dub the 
three models PopIII-Edd, Massive-MH and Massive- 
subEdd respectively. Note that in the model names the first 
part refers to the type of seed and the second part refers to 
the kind of accretion history assumed. Therefore the model 
PopIII-Edd refers to: initial sees from Pop HI remnants 
always accreting at the Eddington rate; model Massive- 
MH refers to initial massive seeds accreting with Edding- 
ton ratios drawn from the MH08 distribution and the model 
Massive-subEdd: initial massive seeds accreting 0.3X Ed- 
dington at all times. 

In the Massive-MH model the accretion rate is not 
limited to the Eddington rate and mildly super-Eddington 
accretion rates (up to /Edd ~ 10) are possible and allowed 
as per MHOS. 

For all three scenarios considered here, accretion starts 
after a dynamical timescale at the virial radius, idyn = 
lO^yr (7?vir/100kpc)(uvir/100kms"^)"\ and lasts until the 
SMBH, of initial mass Min, has accreted AM. The lifetime 
of an AGN therefore depends on how much mass it accretes 
during each episode: 



iAGN 



^Edd 



/Edd 1 






(5) 



where e is the radiative efficiency (e ~ 0.1), 
tedd = 0.45 Gyr and Mfln = min[(Min + AM), 1.3 x 
10* (a/200 kms-i)*-2*M0]. 

We further assume that, when two galaxies hosting 
SMBHs merge, the SMBHs themselves merge within the 
merger timescale of the host halos, which is a plausible as- 
sumptio n for SMBH binar ies formed after gas-rich galaxy 
mergers iDotti et al.l (12007). W e adopt the relations sug- 
gested by iTafFoni et al.l (2003) for the merger timescale. 
Black holes are allowed to accrete during the merging pro- 
cess if the timescale for accretion, corresponding to the sum 
of the dynamical timescale and ^agn, is longer than the 
merger timescale. 

As outlined earlier, in propagating the seeds it is as- 
sumed that accretion episodes and therefore growth spurts 
are triggered only by major mergers. We find that in a 
merger-driven scenario for SMBH growth the most biased 
galajcies at every epoch host the most massive SMBHs that 
are most likely already sitting on the Mbh — cr relation. 
Lower mass SMBHs (below lO'' Mq) are instead off the re- 
lation at 2 = 4 and even at z = 2. These baseline results are 
independent of the seeding mechanism. In the initial 
massive seeds scenario, most of the SMBH seeds start out 
well above the z = Q Mbh — o relation, that is, they are 
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Figure 1. The derived MF of SMBHs at z = 4.25. The up- 
per shaded curve in all three panels is the MF for all SMBHs 
including active and inactive ones. The lower (darker) shaded 
curve in all three panels is the MF for SMBHs that can be 
identified as BLQSOs. The dashed curve in all three panels 
is the MF of all SMBHs from MHOS. The solid curve is the 
MF of BLQSOs from KIO. The three panels refer to the mod- 
els: PopIII — Edd (uppermost panel). Massive — MH (middle 
panel). Massive — subEdd (bottom panel). 



'over massive' compared to the local relation. Seeds form 
only in haloes within a narrow range of velocity dispersion 
[a ~ 15kms~^ at the earliest epochs, see eqns. 1 and 3. The 
SMBH mass corresponding to cr ~ 15kms~^, according to 
the local A/bh — cr relation, would be ~ 3 x 10^ Mq. The 
MF instead peaks at 10^ M© (Lodato & Natarajan 2007). 
As time elapses, all haloes are bound to grow in mass by 
mergers. The lowest mass haloes, though, experience mostly 
minor mergers, that do not trigger accretion episodes, and 
hence do not grow the SMBH. The evolution of these sys- 
tems can be described by a shift towards the right of the 
AfBH — <y relation: a increases, but A/bh stays roughly con- 
stant. 



4 RESULTS AND CONCLUSIONS 

In what follows, we present a detailed comparison of the 
data with our three models, namely the PopIII-Edd: ini- 
tial seeds from Pop HI remnants with accretion assumed at 
all times at the Eddington rate; Massive-MH: wherein the 
initial massive seeds have accretion rates drawn from the dis- 
tribution determined by MHOS; Massive-subEdd: initial 
massive seeds with accretion assumed at all times to be at 
0.3 X the Eddington rate. For all three models, we compare 
our derived MF of BLQSOs with that estimated by KIO and 
the MF of all SMBHs (the entire population that includes 
actively accreting -I- inactive black holes) with the MF pre- 
dicted by MHOS. Although the latter is also a theoretical 
MF, it is derived under completely different assumptions. 
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Figure 2. The derived MF of SMBHs at 2 = 3.25. The up- 
per shaded curve in all three panels is the MF for all SMBHs 
including active and inactive ones. The lower (darker) shaded 
curve in all three panels is the MF for SMBHs that can be 
identified as BLQSOs. The dashed curve in all three panels 
is the MF of all SMBHs from MHOS. The solid curve is the 
MF of BLQSOs from KIO. The three panels refer to the mod- 
els: PopIII — Edd (uppermost panel), Massive — MH (middle 
panel), Massive — subEdd (bottom panel). 



Figure 3. The derived MF of SMBHs at 2 = 2. The upper shaded 
curve in all three panels is the MF for all SMBHs including active 
and inactive ones. The lower (darker) shaded curve in all three 
panels is the MF for SMBHs that can be identified as BLQSOs. 
The dashed curve in all three panels is the MF of all SMBHs from 
MHOS. The solid curve is the MF of BLQSOs from KIO. The three 
panels refer to the models: PopIII — Edd (uppermost panel), 
Massive — MH (middle panel). Massive — subEdd (bottom 
panel). 



as the MHOS model starts with SMBHs today and evolves 
them back in time via a continuity equation. We therefore 
consider this an interesting case of model confronting model 
(see also Volonteri & Begelman 2010). 

We compare our derived MFs for BLQSOs to KIO and 
our MFs for the entire BH population (active -I- inactive) to 
MHOS in Figures 1-4. We find the following: 

• In the higher redshift slices at z — 4.25 and z — 3.25, 
independent of our models the MF derived in MHOS for all 
SMBHs is incompatible with that for the subset of BLQSOs 
from KIO. In these instances the abundance of BLQSOs from 
KIO is in excess of that of all SMBHs from MHOS models. 
It is only in the lowest redshift bin at 2: = 1.25 that the 
MFs of BLQSOs from KIO is lower than that of all SMBHs 
derived from the MHOS models. It is therefore clear that the 
Eddington ratio distributions derived from KIO and from 
MHOS are most likely discrepant. 

• At the highest redshift slice z — 4.25, the light seed 
model PopIII-Edd) shown in the top panel fails to grow 
SMBHs massive enough (more massive than ~ 10^'^Mq) 
to populate the high-mass end of the MF of BLQSOs of 
KIO. The discrepancy is roughly two orders of magnitude as 
clearly seen in the upper panel of Figure 1. Since the Ed- 
dington ratio was assumed to be unity to evolve this model, 
a deficit at the high mass end cannot be compensated for 
by altering the accretion physics any further. Alternatively, 
if SMBHs can become "over massive" at higher z (i.e., if 
they climb well above the standard Mbh — o" relation) and 
have phases of super-Eddington accretion (Volonteri & Rees 



2005), then enough SMBHs at the high mass end might be 
produced. We note that assuming the MHOS accretion model 
in conjunction with Pop HI remnants (this model is not 
plotted in our figures) makes the discrepancy at the high 
mass end much worse, as the accretion rate for low mass 
black holes (< 1000 Mq) in the MHOS model peaks at sub- 
Eddington rates at all redshifts. 

• Both the massive seed models (Massive-subEdd and 
Massive-MH) also under-produce the KIO MFs at the 
highest redshift slice as seen clearly by looking at the middle 
and bottom panels in z = 4.25; 

• The MF of all existing SMBHs (active + inactive) 
from the PopIII-Edd, Massive-subEdd and Massive- 
MH models overestimates MHOS's MF for all SMBHs at all 
masses > 3 x 10^ A/0 at z = 3.25, while there is reason- 
ably good agreement with KIO for the BLQSO MF for the 
Massive-MH model. 

• At lower redshifts {z — 1.25 and z = 2) both massive 
seed models start to overestimate the MF of BLQSOs by 
KIO at Mbh > 10^ Mq. This is not unexpected, as we have 
not implemented a cut-oflto mimic the depletion of available 
gas in massive galaxies. 

• All massive seeds models (Massive-subEdd and 
Massive-MH) overestimate the MF of all SMBHs (active -|- 
inactive) from MHOS at Mbh > 10*'^ Mq. At 2 < 2 we can- 
not clearly assess whether MHOS might underestimate the 
total MF. We are, however, inclined to attribute the discrep- 
ancy to our models, for the same reason that we likely over- 
estimate the MF of quasars as we do not take into account 
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Figure 4. The derived MF of SMBHs at 2 = 1.25. The up- 
per shaded curve in all three panels is the MF for all SMBHs 
including active and inactive ones. The lower (darker) shaded 
curve in all three panels is the MF for SMBHs that can be 
identified as BLQSOs. The dashed curve in all three panels 
is the MF of all SMBHs from MHOS. The solid curve is the 
MF of BLQSOs from KIO. The three panels refer to the mod- 
els: PopIII — Edd (uppermost panel), Massive — MH (middle 
panel), Massive — subEdd (bottom panel). 



changes in the gas inventory at late times in the universe. 
The PopIII-Edd model gets instead into better ahgnment 
at lower and lower redshifts. At z = 2 these MFs begin to 
agree up to 10* Mq and at 2 = 1.25 they fall into excellent 
agreement up to 10^ Mq . 

• For model Massive-MH we can estimate the typical 
Eddington ratio of SMBHs powering BLQSOs. We find that 
most luminous quasars at the highest redshift have accretion 
rates close to or even slightly larger than the Eddington rate 
(the fraction of super-Eddington accreting SMBHs is ~ 60% 
at « = 3; ~ 40% at 2 = 2; ~ 20% at 2 = 1). At lower 
redshift, the flux limit of the SDSS survey corresponds to 
lower luminosities, and less powerful accretors are selected. 
In general, we find that the accretion rates in the models 
overestimate the accretion rates derived by KIO. 

• We estimate the duty cycle, defined as the fraction of 
SMBHs that are active and detectable as BLQSOs, and show 
the results for the three models (PopIII-Edd, Massive- 
subEdd, Massive-MH) in Figure [5] at 2 = 1, 2 = 2, 2; = 
3, 2 = 4 (dark to hght shade of colours). Comparing the 
theoretical duty cycle with the duty cycle estimated by KIO 
at 2 = 1, we find that model PopIII-Edd is consistent at 
all masses; while models Massive-MH and Massive-Edd 
overestimate the duty cycle at A^bh > 10^ Mq (as expected 
from the MF shown in Figures 3 and 4.). 

• We find that the duty cycle increases with increasing 
redshift for all our models, although more strongly for the 
Massive-subEdd model. This is due to the fact that this 
model is the one with the lowest average accretion rate 



(/Edd = 0.3, versus /Edd ~ 1 for the PopIII-Edd and 
/Edd^l in Massive-MH at 2 > 3, therefore, from Equa- 
tion [5] the BHs in the Massive-subEdd model need to ac- 
crete for a longer time in order to reach the final mass set 
by the Mbh — cr relation. 

In conclusion, the implications for seeding models are 
that light seeds always under-produce the bright end of the 
BLQSO LF for high black hole masses. This is a serious issue 
as they are assumed to be already accreting at the optimal 
Eddington rate. 

Moreover, according to recent work (Turk et al. 2009; 
Greif et al. 2011; Davis et al. 2011) it appears now that (i) 
the IMF of the first stars may not be biased high as frag- 
mentation might occur during the formation process and (ii) 
due to turbulence the formation of the first stars might get 
delayed. So the PopIII-Edd model aside from being unable 
to match the black hole MF derived from BLQSOs from the 
SDSS as shown here, might in fact not be the most efficient 
channel to produce seeds. It is clear that a channel to form 
more massive black hole seeds at high redshifts is required 
to start matching observations of the bright end (massive 
end of the black hole MF) from the earliest to intermediate 
redshifts (2 ~ 1), and not only at the highest redshifts as 
previously noted (e.g., Haiman 2004; Yoo & Miralda-Escude 
2004; Shapiro 2005; Volonteri & Rees 2005, 2006). There- 
fore, low-mass SMBH seeds such as Pop III star 
remnants require exceptional growth conditions not 
only to explain the existence of rare objects such as 
2 > 6 quasars, but also to explain the unexceptional 
population of luminous quasars at 2 ~ 2 — 3. 

Over-predicting the BHMF at the high mass end is typ- 
ically not a serious problem as lowering the Eddington rate 
would fix that. For instance, Natarajan & Treister (2009) 
propose the existence of an upper mass limit to BHs at all 
epochs, arising from the shutting off of accretion due to self- 
regulation of the gas supply in the galactic nuclei. According 
to their estimate this upper limit is given by: 



Mbh ~ 5 X 10^ (- "" 



350 km s 



'Ma 



(6) 



KIO estimate an upper mass limit at 2 ~ 1 for a black hole 
hosting a BLQSO to be Mbh ~ 3 x 10^" Mq. Utilizing the 
limit proposed by Natarajan & Treister (2009), if a BH shuts 
off at the end of this BLQSO episode at 2 ~ 1, it is likely 
hosted in a galaxy with a velocity dispersion of a ~ 440 
km/s, which corresponds to the typical velocity dispersion 
of the brightest central galaxy in a cluster. This fits in nicely 
with the view of down-sizing wherein accretion is noted to 
shut off at 2 ~ 1 in these BCGs in clusters. 
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Figure 5. The duty cycle of BLQSOs as a function of SMBH 
mass. The models shown here are again in the same order as be- 
fore: PopIII — Edd (uppermost panel), Massive — MH (mid- 
dle panel), Massive — subEdd (bottom panel). In each panel, 
we show the duty cycle at z = 4, z = 3, z = 2, and z = 1 (from 
lightest to darkest shade respectively). The square symbols show 
the lower bound on the duty cycle at z = 1 from KIO. 
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